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ABSTRACT 

We study the expected properties and statistics of giant arcs produced by galaxy clusters in a ACDM 
universe and investigate how the characteristics of CDM clusters determine the properties of the arcs 
they generate. Due to the triaxiality and substructure of CDM halos, the giant arc cross section 
for individual clusters varies by more than an order of magnitude as a function of viewing angle. 
In addition, the shallow density cusps and triaxiality of CDM clusters cause systematic alignments 
of giant arcs which should be testable with larger samples from forthcoming lensing surveys. We 
compute the predicted statistics of giant arcs for the ACDM model and compare to results from 
previous surveys. The predicted arc statistics are in excellent agreement with the numbers of giant 
arcs observed around low redshift (0.2 < z < 0.6) clusters from the EMSS sample, however there are 
hints of a possible excess of arcs observed around high redshift z > 0.6 clusters. This excess, if real, 
appears to be due to the presence of highly massive or concentrated clusters at high redshifts. 
Subject headings: gravitational lensing — galaxies: clusters — dark matter 



1. INTRODUCTION 

Clusters of galaxies provide some of the most spec- 
tacular examples of gravitational lensing. The popu- 
lar images of Abell 2218 and CL0024+1654 that are 
prominent on the Hubble Space Telescope public website 
(www.hubblesite.org) show very elongated multiply im- 
aged background galaxies whose features have been dis- 
torted by the deep gravitational potential of the galaxy 
cluster. Such cluster lenses have a variety of cosmolog- 
ical uses. For example, these natural gravitational tele- 
scopes can greatly magnify distant sources, allowing us to 
stud y the proper ties of otherwise unobserv able galaxies 
Ce.g. IBlain et al.lll99a iMetcalfe et 311120031: ISmail et all 
l2002j) . In addition, strongly lensed arcs offer a unique, di- 
rect probe of the cluster gravitational potential on scales 
where dark matter is expected to be the dominant com- 
ponent. Heating, cooling, and star formation can have a 
significant impact on baryons but the dark matter, and 
therefore the total potential, should be relatively unaf- 
fected by these poorly understood processes. Lastly, the 
incidence of giant arcs may be used to s tudy the back- 
ground cosmology itself. iBartelm ann et al. (1998, here- 
after B98) found that the predicted number of giant arcs 
varies by orders of magnitude among different cosmolog- 
ical models. In light of mounting evidence supporting 
the ACDM model (e.g. Spergel et al. 2003), the spectac- 
ular images of giant arcs have led to an embarrassment 
of riches, as B98 found that the observed instances of 
giant arcs exceeded their predicted rate by an order of 
magnitude. 

The original discrepancy reported in B98 was based 
on a subsample of 16 of the most massive clusters 
l|Le Fevre et al.lll994D in the EMSS survey and the arc 
frequency (roughly 20% of very massive clusters) was 



confirmed in the larger EMSS sample of iLuppino et alJ 
(1999). A subsequent survey based on optically-selected 
clusters in the Las Campanas Distant Cluster Survey 
(jZaritskv fc Gonzalez 2003) found comparable giant arc 
frequencies, while a recent re port by the Red Cluster Se- 
quence cluster survey group ijGladders et 31.1120031 12002D 
found a high probability (~ 30%) of a lensing cluster 
showing more than one source being distorted into a gi- 
ant arc. 

On the theoretical side, significant work has gone 
into refining the expected number of giant arcs. Sev- 
eral works have confirmed th e lensing cross-sections of 
B98 jBartelmann et all 120031: iMeneghetti et alJ l2003bl 
120011) . but other recent work has suggested that the 
cr oss-sections of B98 may be too low. In particu- 
lar [^^^tegans^^Q l|2003(l used numerical simulations 
and found that the cross-section was a very steep func- 
tion of redshift and that the anomalous cross-sections 
could be brought into agr eement by allowing a broader 
range of source redshifts. lOguri et al . (2003) used ana- 
lytic models with triaxiality and found that allowing a 
steeper central density profile enhanced the lensing cross- 
section to a level that was close to the observed lensing 
frequency of the EMSS sample. 

In this paper we go back to the beginning of the prob- 
lem (i.e., B98) to repeat the analysis on a larger sam- 
ple of simulated clusters and update the statistics for 
the EMSS lensing sample. We then apply the analy- 
sis to recent optical catalogs to determine the extent, if 
any, of the discrepancy between observations and the- 
ory. Lastly, we investigate what new information may 
be gleaned from upcoming deep, large area surveys like 
the CFHT Legacy survey, which should produce large 
numbers of new strongly lensed arcs. 
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2. CLUSTER LENSING SIMULATIONS 
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We estimate giant arc statistics by ray-tracing through 
cosmological N-body simulations. Below, details of the 
calculation are briefly described. 

2.1. N-body simulations 

We used the publicly available simulation outputs 
of the Virgo Project^. We used the A CDM simula- 
tions that were performed as part of the GIF project 
(pCauffmann et al. 1999). The cosmological parameters 
were il™ = 0.3, JIa = 0.7, h = 0.7, power spectrum shape 
parameter T = 0.21, and mass power spectrum normal- 
ization fjg = 0.9. The comoving simulation box length 
was 141.3 Mpc and 256'^ particles were used, leading 
to a mass per particle of 10^^h~^MQ. The gravitational 
softening length was 20/i~^kpc. Simulation outputs were 
available at a variety of redshifts, allowing accurate al- 
lowance of time evolution in cluster properties. A subset 
of the clusters from these simulations were used as part 
of the ensemble of simulated clusters used in B98. 

2.2. Lensing calculation 

To compute giant arc statistics for the N-body sim- 
ulations described above, we adopt a technique similar 
to that pioneered by bartelmann fc Weiss ( 1 99^ and 
empl o yed by subsequent worke r s (e.g Bartel mann et al.l 
IT9981 IMenes-hetti etaP I2003ai I2000L 2001.') . We first 
project the dark matter density onto a two dimensional 
regular grid. Typically we used a 256 x 256 grid covering 
a square Ah~^ comoving Mpc on a side, centered on the 
cluster. From the projected surface density E, the con- 
vergence for each lens and source redshift pair is given 
by K = Yi/Yjcrit{zh Zs), wh ere Sf-rit. is the critical surface 
density for strong lensing (I Schneider et al.lll992(l . Then, 
the deflection angle map cy.{x) is constructed (via Fourier 
transform) on the same grid. Given this map, we then 
ray trace a finer 4000 x 4000 grid covering a subset of the 
image plane onto the source plane, linearly interpolating 
the deflection between the grid points of the coarser 256^ 
grid. We set up a regular grid on the source plane; each 
source plane pixel contains a linked list of all image plane 
pixels which ray trace to the source plane pixel. Having 
performed this single ray trace, it is then a simple mat- 
ter to construct the lensed arcs for an arbitrary source by 
concatenating the lists for all source plane pixels falling 
in the source. We emphasize that the use of linked lists 
greatly speeds the calculation ~ only one single ray trace 
is performed, after which arcs may be constructed for 
many different source positions, sizes and geometries for 
essentially no computational cost. 

For each source, the pixel lists are then sorted into 
separate arcs, by grouping together neighboring pix- 
els. We measure arc proper ties using methods based on 
iBartelmann & WeissI l(T99^ . The arc area and magni- 
fication are found by summing the areas of the pixels 
falling in the arc. Arc lengths are estimated by first find- 
ing the arc centers, then finding the arc pixel furthest 
from the centroid, as well as the pixel furthest from this 
pixel. As in B98, the arc length was then given by the 
sum of the lengths of the two line segments connecting 
these three points. The arc width was defined as the ra- 
tio of the arc area to the arc length. See figure Q for an 
example. 

^ http:/ /www.mpa-garching.mpg.dc/Virgo 




Fig. 1. — Example of simulated arcs. The source plane (left) 
and image plane (right) are illustrated. Solid lines depict the ap- 
proximate locations of the critical lines and caustics, while filled 
regions show the pixels falling inside the source (left panel) and the 
corresponding lensed images (right panel). Two arcs are found for 
this source, with length to width ratios of 7.5 and 2.8 respectively. 
For the long arc, the arc radius r and position angle (relative to 
the major axis of the mass distribution) are labeled. 
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Fig. 2. — Magnification vs. length/width ratio. While these 
two quantities are equal for isothermal lenses, giant arcs tend to 
form at small radii (r < Ts) where the density profile is shallower 
than isothermal. Here we plot all arcs in one cluster projection 
with magnification /i > 5, however only ~ 40% of the points have 
length/width ratios r > 5. 



The relation between length/ width ratio and magni- 
ficati on is sensitive to t he slope of the halo potential 
( Willi ams fc Lewi^ll998j) and can also be affected by lo- 
cal perturbations in the projected mass distribution. In 
figure 121 we plot arc magnification against length/ width 
ratio using the above method. We find good agreement 
for the mean trend o f length/width vs. magn ification 
with the predictions of lWilliams Sz Lewi^ l|1998() for typ- 
ical clusters in the simulations. However, there is sub- 
stantial scatter in this relation that is likely caused by 
random source orientations and local fiuctuations in the 
surface mass density. It is apparent that highly distorted 
images will also be highly magnified, but the converse is 
not true. 

For each cluster projection, we compute giant arc 
statistics using a Monte Carlo approach. We randomly 
place elliptical sources, with axis ratios drawn uniformly 
in the range q £ [0.5, 1] and random position angles. We 
count as giant arcs only those images with length/ width 
ratios r > 10, making no distinction between radial and 
tangential arcs, although we note that the latter dom- 
inate the cross section. Cross sections are then com- 
puted by multiplying the fraction of sources producing 
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Fig. 3. — Histogram of cross-sections for the 4*^ most massive 
cluster at 2 = 0.2 under different projections. 



giant arcs in the Monte Carlo by the area over which the 
sources were randomly distributed. 

2.3. Cross sections 

The giant arc cross section for a single cluster varies 
enormously as a function of orientation, as shown in 
fig- 01 A factor of ~ 20 spread in a is not uncommon. 
Additionally, the distribution is quite skewed, with a tail 
extending to high a. This appears to be due to the large 
degree of substructure and ellipticity/triaxiality in the 
density profiles. Orientations with large cross sections 
are typically those which project along the cluster major 
axis, or which project massive substructures onto small 
radii. Unfortunately, this means that shot noise will be 
persistent in any calculation of arc statistics; large num- 
bers of projections are required for trustworthy mean 
cross sections. In the calculations below, we compute 
mean cross sections by averaging over 251 orientations 
for each of the top 20 most massive clusters at each red- 
shift, while we perform 3 projections each for the next 
100 clusters. 

In addition to the scatter in cross section for individ- 
ual clusters, there is large scatter among cross sections 
of clusters of similar mass, shown in Figure 0| The error 
bar for each cluster shows the error in the mean, averaged 
over all orientations. The dispersion for each cluster is 
comparable to the mean. The mean lensing cross-section 
as a function of mass shows a dramatic rise at a mass 
of a few lO^** h~^MQ and an apparent flattening at the 
highest masses; the numbers become too low to say any- 
thing definite at high masses, as shown by the large error 
bar in the highest mass bin. 

When weighted by the number of clusters at each 
mass (right panel) there is an apparent peak at ~ 7 x 
10^'^h~^MQ. The significance of this peak is hard to 
judge; given the low numbers of clusters in the highest 
mass bin, and the large scatter among cross sections in 
each bin, we cannot definitively say whether we are un- 
derestimating the total optical depth by neglecting mas- 
sive clusters too rare to appear in our limited simulation 
volume. On the other hand, the peak corresponds ex- 



actly to the mass range expected of massive X-ray lu- 
minous galaxy clusters which appear to dominate the 
giant arc optical depth observed by surveys like EMSS. 
In addition, we know that the rise of (a) with m must 
flatten at high masses, once the Einstein radii approach 
the NEW scale radius Vg. This ensures that, at some high 
mass, the exponentially declining mass function n(rn) 
will overpower the growing {a)(m). We will assume that 
our simulation volume provides a fair sample of the ob- 
jects causing most of the lensing, with the caveat that 
our estimates might plausibly be low by a factor of ~ 2. 

A related caveat is that, since the total optical depth 
is dominated by the most massive clusters, great cau- 
tion must be employed before deriving cosmological con- 
straints from giant arc statistics. While the abundance 
of such clusters depends in part on cosmological parame- 
ters like r^A, it also depends sensitively upon the matter 
powe r spectrum fe.g. lBahcall &: FanlllQQStlHaiman et alJ 
I2001D . Strong prior constraints on parameters like as 
must be applied (and believed) before constraints can be 
derived on parameters like Q\ or w. 

2.4. Arc properties and cluster galaxies 

The effects of cluster galaxies on the formation of giant 
arcs have been investigated bv iMeneehetti et all f2000) 
and lElores et al. (2000.} . and the effe cts of central cD 
galaxie s in particular were studied bv iMeneghetti et all 
()2003a|) . These authors found that cluster galaxies are 
generally unimportant, except for central cD galaxies 
which can enhance the giant arc cross section by 50%. 

The relatively small effect of cluster galaxies on the 
arc formation cross section is not surprising. Giant arcs 
typically form at large radii relative to the cluster cen- 
ter, ~ 20 — 30"; see figure A mass m « Scrit • « 
lQ^'^h~^MQ is interior to these arcs, much larger than the 
typical masses even of cD galaxies, < 10^^ Mq. Accord- 
ingly, we do not expect the cross sections for formation 
of these wide separation arcs to be grossly affected by 
the presence or absence of galaxies in the simulations. 
On the other hand, galaxies can strongly affect the cross 
section for formation of smaller separation arcs. To il- 
lustrate, we plot in figure IHl the distribution of arc radii 
produced by clusters from the z = 0.3 output, lensing 
sources at = 2. We artificially added mass concentra- 
tions to the cluster centers in two ways: as point masses 
of varying mass, and as galaxies with fixed velocity dis- 
persion (350 km/s) and varying mass. It is apparent from 
the figure that mass is a less important determinant of 
arc properties than central mass concentration, here pa- 
rameterized by velocity dispersion. Galaxies of different 
masses (e.g. 3 x 10^^ and IO^'^/i^^M^q) but the same ve- 
locity dispersion produce nearly indistinguishable effects 
on the radial distribution. Increasing the central mass 
concentration can greatly increase the number of giant 
arcs forming at small radii, illustrated by the curves for 
which point masses were added to the central pixel. How- 
ever, it would appear that arc statistics at very large radii 
(> 20") are unaffected by these central mass concentra- 
tions. Even at moderate radii (> 10"), the difference 
between the statistics for pure dark matter simulations 
and those including reasonably sized galaxies is not dra- 
matic. 

Accordingly, we expect that our pure dark matter sim- 
ulations should provide a reasonably accurate estimate 
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Fig. 4. — Cross-section for forming giant arcs (loft) and cross-section weighted by number of clusters (right) for simulation volume at 
zi = 0.2 and a source redshift Zs = 1. Cross sections are given in physical distances on the source plane. The solid line in the left panel is 
the mean cross section (cr) averaged over bins of width 0.2 in logj^Q M. Error bars in the left panel show error in the mean cross-section 
per cluster, not the dispersion, while the error bars in the right panel include the effect of Poisson statistics. 
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Fig. 5. — Radial distribution of giant arcs, as measured by 
ILuppino et a l. (1999) (solid line), and as predicted by the z = 0.3 
GIF clusters (dashed line). 
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Fig. 6. — Number of Zs = 2 giant arcs produced by a (141 
/i~^Mpc)^ volume of clusters at z = 0.3. The solid line shows 
results from the dark matter distribution alone, the 'P' curves show 
results when a central point mass is added, and the 'C curves 
show the effects of adding a central singular isothermal sphere with 
velocity dispersion a = 350 km/s. 



for the optical depth for wide separation giant arcs in 
the ACDM modeL However, absent a prescription for 
galaxy formation, we cannot make reliable predictions 
for the numbers of intermediate separation (r < 5") arcs, 
which will prove to be a limitation in the interpretation 
of high redshift giant arc surveys. 

Although cD galaxies do not dramatically change the 
wide separation arc cross section, they can have im- 
portant effects on ^ant arc properties. For example, 
l^illiams ct al. (1923) have argued that massive (M ~ 
3 X IO^^Mq) galaxies are required at the centers of most 
low redshift arc-bearing clusters in order to explain ob- 
served arc widths. Central galaxies also affect the an- 
gular distribution of long, thin arcs about the cluster 
center. As an example, we plot in figure |7| the angular 
distribution of giant arcs about one particular cluster, as 
a function of the mass of the central cD galaxy painted 
onto the surface density. In the top panel, corresponding 
to pure DM with no added galaxy, the giant arcs tend to 
form only at the ends of the critical line. As the mass be- 
comes more centrally concentrated, the critical lines do 
not change dramatically in size, however giant arcs begin 
to form all around the critical lines, instead of only along 
the cluster major axis. 

This is shown more quantitatively in figure |S| which 
plots the angular distribution of arcs about the cluster 
center, relative to the cluster major axis. About 80% 
of giant arcs form at position angles within 45° of the 
cluster major axis in the pure DM simulations. The 
large anisotropy in the angular distribution is diminished 
somewhat by the addition of the central galaxy; never- 
theless we expect a residual anisotropy to be present for 
realistic galaxy masses. For example, adding a central 
mass of 3 X 10^^ Mq decreases the fraction of arcs 
within 45° of the major axis to ^ 70%. Under the as- 
sumption that the cluster major axis can be determined 
from the brightest central galaxy position angle, this pre- 
diction can be tested using clusters even which have only 
one arc. This anisotropy signal may be washed out by 
misalignments between the BCG orientation and halo 
orientation, however for multi-arc clusters the anisotropy 
should also be manifest in the distribution of relative po- 
sition angles of arcs, 4>i2 = 4'i — 4>2, also shown in figure|H| 
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Fig. 7. — Angular distribution of giant arcs for one cluster projec- 
tion. The panels correspond, from top to bottom, to results from 
the DM simulation with added central point masses of 0,1, and 3 
10^^h~^ Mq respectively. The solid lines are the critical curves, 
and symbols show the centers of giant arcs produced in the Monte 
Carlo. 
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Fig. 8. — Distribution of arc position angles from z = 0.3 GIF 
clusters. The left panel shows arc position angles relative to the 
major axis of the projected dark matter distribution, for pure DM 
(solid) and DM + 3 X 10^^ h~^MQ (dashed). The histogram shows 
the same distribution for the EMSS arcs. The right panel shows 
the distribution of relative position angles between pairs of arcs in 
multi-arc systems, 0i2 = (^i ~ <t>2, for DM only. The histogram 
shows the </>i2 distribution for the three EMSS clusters showing 
multiple giant arcs. 



This behavior can be easily understood, as outhned 
in more detail in the appendix. The caustic structures 
of the pure dark matter case often take the form of 
"lips" catastrophes, while the addition of the central 
mass transforms these into the more familiar radial and 
tangential caustics. Lips caustics tend to form giant arcs 
only at the ends of the critical line, while tangential caus- 
tics can form giant arcs at a wide variety of position an- 
gles. 

Since we have no reliable model for populating CDM 
halos with galaxies, we cannot at present make predic- 
tions regarding the angular distribution of arcs around 
clusters. However, our crude treatment of central galax- 
ies seems to indicate that typical cD galaxies do not fully 
isotropize the arc angular distribution. It will therefore 
be interesting to measure the angular distribution with 
upcoming wide-angle, deep lensing surveys which should 
find numerous multi-arc clusters. If clusters have lips. 
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Fig. 9. — Histogram of cross-sections for the same M = 
4 X 10^'^h~^ Mq cluster at 2 = 0.3 simulated at two different res- 
olutions. The shaded histogram corresponds to the high resolu- 
tion simulation, with a mean cross section for producing l/w > 10 
arcs at d > 10" of (a) = 210(/i-lkpc)2. The unshaded his- 
togram corresponds to the low resolution simulation, which gives 
(ct) = 250(/i-^kpc)2. 



then multiple arcs in a cluster should tend to form along 
the cluster major axis, at the ends of the critical curve. 

2.5. Numerical resolution 

We argued above that our dissipationless N-body sim- 
ulations cannot be used to make reliable predictions for 
arc statistics at small radii, while at large radii {d > 
10 — 20"), neglect of baryonic physics should not change 
our results significantly. Even at large radii, however, an- 
other possible source of error is particle noise in the simu- 
lation. Collapsed objects represented by low numbers of 
particles will suffer artificially from two-body relaxation 
effects, which can si gnificantly aff e ct the inner structures 
of simulated halos. iPower et al.l l)2003j) find that halos 
satisfying Wx(»^) = idyniV(< r)/[81ogiV] > O-GHq-^ 
may safely be assumed to be free of relaxation effects. 
We are interested in measuring arc properties at radii r ~ 
50 - lOO/i-^kpc for objects of ma ss M > IO^-^/i-^Mq. 
For such objects, the IPower et al.| bound implies a mini- 
mum of A^(rEin) ^ 1500 particles representing the profile 
interior to the arc radius, or a particle mass not exceeding 
mp < 4 X IO^/j-^Mq. For comparison, the ACDM GIF 
simulations have particle masses of 1.4 x 10^^h~^ Mq, in 
excess of this bound. Hence, it is unclear from this simple 
analysis if our results will be affected by particle noise. 



Two-body relaxation can also become imp ortant if the 
force softening is insufficient. IPower et al.l argue that a 
force softening length e « 4rvir/-\/-^vir is optimal. For our 
example halo this works out to about eopt ~ 16/i~^kpc; 
while in comparison the GIF ACDM simulations use a 
softening of 20h~^ kpc. 

To test whether such resolution effects will signifi- 
cantly impact on our conclusions, we perform the follow- 
ing convergence test. We perform two simulations of a 
{100h''^Mpc)^ volume at two different resolutions, using 
the TreePM code^ (kindly made publicly available by P. 
Bode and J. P. Ostrikcr). The low resolution simulation 

see ]http://www.astro.princeton.edu/~bode/TPM/index.html| 
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uses 128'^ particles with mass irip = 3.6 x lO^^h ^Mq 
and force softening length e = 2Qh^^ kpc, while the 
high resolution simulation has 256^ particles of mass 
rup = 4.5 X IO^H^^Mq and softening e = 10h~^ kpc. 
To ensure that the same objects would form in both sim- 
ulations, they were initialized using the same random 
phases for all modes common to the two boxes, while the 
higher resolution run contained additional high-fc power. 
We then compared individual clusters from each simu- 
lation. One example is shown in Figure the other 
clusters we checked showed similar levels of agreement 
between the two resolutions. We therefore conclude that 
the resolution of the GIF simulations is adequate for our 
purposes. 

3. HOW MANY LUMINOUS GIANT ARCS? 

Following the considerations of ij2 2.41 we will count in 
our statistics only those arcs with length-to- width ratios 
r > 10, and distances relative to the cluster center of 
d > 10" in order to minimize the uncertainties from the 
galaxy contribution. 

The optical depth for forming an arc can be derived 
for each simulation volume as a sum over all clusters in 
the box: 

dr _ dVco/dzi ^jjzs) . . 

The total optical depth as a function of source redshift 
t{zs) is then simply a sum over simulation volumes. The 
number of arcs expected per square degree is then given 

by 

f°° dNsource 
Narcs = / dZsT{Zs) Epill ^ (2) 

Jo dzs 

where dNsource/ dzs is the source density that could form 
observable arcs. 

In figure [TUI we plot the contribution to the giant arc 
optical depth from different lens rcdshifts, for three dif- 
ferent source redshifts. As Wambsganss et al. (2003) 
have stressed, the optical depth increases with source 
galaxy redshift: we find optical depths of 0.25, 0.70, & 
1.4 X 10~® for source redshifts Zg — 1, 1.5, and 2 respec- 
tively. However, we find that the variation with redshift, 
while substantia l, is not as strong as that reported by 
iWambseanss et all (|2063i) . 

The source density term is of equal importance to the 
optical depth term. To estimate the number density of 
galaxies as a function of redshift and limiting magnitude, 
we em ploy the photometric catalogs of Fontana et al] 
l)200n() . compiled from galaxy counts from the Hubble 
Deep Fields. We divided source galaxies into three bins: 
0.75 < < 1.25, 1.25 < < 1.75, and Zs > 1.75. 

Magnification bias may or may not be important. In 
a strictly flux-limited survey, magnification bias would 
be important, but there is no magnification bias in the 
limit where all sources are resolved and surface bright- 
ness is the limiting concern. For giant arc surveys, nei- 
ther regime is appropriate. Arc detection is closer to 
being limited by surface brightness than by integrated 
fiux, but adaptive smoothing can modify this. The ob- 
served source distribution is expressed in terms of inte- 
grated fiux, but ground-based surveys may be missing 



Table 1. Survey statistics 



Survey 


redshifts 


area 


depth 




-^prcd 


EMSSi 


0.15<^<0.6 ' 


- 360 


V <22 


8 


8.2 


LCDCS2 


0.5 <z < 0.7 


69 


R < 21.5 


2 


1.2 


RCS^ 


z < 0.6 


90 


Mfl < 24 





2 




z > 0.6 






4 


1 


Note. — 


Survey areas are 


given 


in square 


degrees. 


For the 



LCDCS, quoted arc statistics are for all radii, with no cut on 
arc radius, while RCS statistics are for arcs in their pr i mary sam- 
ple with d > 10". Referen ces a re: fl) [Li^ppino etaljl 119931 . (2) 
IZaritskv &: Gonzaled <2003) . CS) IGladders et al.l I200'3i . 



distant compact sources because seeing can reduce their 
apparent surface brightness. In what follows we include 
magnification bias when comparing to catalogs that are 
nominally magnitude-limited. 

3.1. The EMSS Sample 

In comparing to observations it is essential that the 
effective area of the survey is well-known. For the 
EMSS, calculating the effective survey area is non-trivial. 
The survey area is a strong function of limiting fiux 
iHenr v etab 1992) . For each cluster detected in the 
EMSS we found the sky area Ai over which that cluster 
could have been detected. An estimate of the cluster sur- 
face density is then given by ^ 1/Ai clus ters per square 
degree. The EMSS lensing sample of 'Luppin o et alJ 
L1399 ) is defined by clust ers with > 2 x lO'^'^erg s~-^ 
and z > 0.15. Using the IHenrv et a l. (1992) coverage, 
we find that the surface density of EMSS clusters is one 
cluster per 9.6 square degrees, a total coverage for the 38 
clusters of 360 square degrees. The earlier sub-sample of 
16 clusters, used in Le Fevre et al . (1994), is defined by 
Lx > 4 X lO^^erg s""'^ and z > 0.2, for which the surface 
density is one cluster per 20 square degrees. 

Of the 38 most massive EMSS clusters, 8 show giant lu- 
minous arcs. Extrapolating over the entire sky we would 
therefore expec t roughl y 900 giant arcs. The sub-sample 
used bv lLe Fevre et al.l 1)19941) showed 6 clusters with gi- 
ant arcs, which translates into one arc-inducing cluster 
every 53 square degrees or roughly 800 giant arcs on the 
entire sky. The consistency of these two values supports 
our estimate for the area covered by EMSS. 

We can roughly estimate the expected number of arcs 
by taking the optical depths from figure 111)1 and multiply- 
ing by an appropriate source galaxy density. Assuming 
that arcs with length/width ratios r > 10 are magnified 
by abo ut 2.5 magnitudes, and noting that Lc Fevre et aO 
(1994) counted arcs with integrated magnitudes V < 22, 
we find galaxy densities for V < 24.5 of 7.4, 1.4, and 
2.6 arcmin^^ for our three redshift bins from the Hubble 
Deep Fields. Multiplying by the corresponding optical 
depths (0.25, 0.70, & 1.4 x 10"^) gives a total of - 940 
arcs expected over the full sky. An alternative estimate 
for the number of arcs may be obtained not by assum- 
ing arcs are magnified by ~ 2.5 magnitudes, but instead 
convolving the calculated arc magnification distribution 
with the observed luminosity function. This estimate 
would be valid in the limit that low surface brightness 
arcs with total integrated magnitudes brighter than the 
detection threshold are observable. For a threshold in- 
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Fig. 10. — The left panel shows the optical depth for arcs with length/width > 10 and distance from cluster center d > 10", as a function 
of redshift for GIF simulations. The right panel shows the number of arcs with total integrated magnitude y < 22 as a function of lens 
redshift, for three different source redshift bins. The declining number density of sources with increasing redshift somewhat counteracts 
the increase of optical depth with redshift. 



tegrated m agnitude of V < 22, t he HDF galaxy counts 
reported bv lFontana "eTal] (l2nonl) give ~ 1400 giant arcs 
with r > 10 located d > 10" from the cluster center, over 
the full sky. 

Apparently, there is no excess in the number of giant 
arcs observed around EMSS clusters compared to expec- 
tations from the ACDM model. It appears that the so- 
lution to what was thought to be an order of magnitude 
problem is in a conspiracy of small effects: a larger effec- 
tive area for the EMSS than estimated in B98 (factor of 
2), a higher source density (factor of 2), and a somewhat 
higher redshift-averaged optical depth (factor of 2). 

3.2. LCDCS Sample 

iZaritskv fc Gonzaled 1)200.'^ surveyed clusters in the 
range 0.5 < z < 0.7 over a 69 deg^ field. In this area, 
they found two giant arcs with R < 21.5 and length 
L > 10", and a number of additional arcs and arc candi- 
dates with which they do not include in their statistics. 
A straightforward comparison with our expected statis- 
tics is complicated by the fact that both of these arcs 
are located close (< 10") to cluster galaxies. As men- 
tioned in W2 2A\ the optical depth for such close-in arcs 
is strongly affected by the presence of galaxies which are 
neglected in the pure dark matter simulations we have 
used. 

Nevertheless, we can try to make some rough estimates 
for the expected giant arc rate. Dropping the d > 10" 
cut we imposed for the EMSS sample, we find that over 
the range 0.5 < z < 0.7, the GIF clusters produce 0.7 
arcs in the 69 deg^ surveyed in the LCDCS. If we paint 
galaxies with a = 350 km/s and A/ = 3 x lO^^/i'^M© 
onto the clusters, the expected number increases to 1.2 
giant arcs. 

Accordingly, it appears that the giant arc incidence 
measured by the LCDCS is consis tent with predictions 
from the ACDM model. However, IZaritskv fc Conzaled 
emphasize that their measured giant arc statistics should 
be viewed as a lower bound on the true incidence; a less 
conservative cut on arcs in their sample could lead to a 
statistically significan t discrepancy. On the other hand, 
IZaritskv &: Gonzaled also note that their estimated clus- 



ter number density is 4.5 x higher than that of the EMSS 
at these redshifts, which could artificially enhance the 
LCDCS arc rate relative to the full sky average. 

3.3. RCS Sample 

IGladders et al] survey 90 deg^ for giant arcs, 

counting only arcs with surface brightness in Rc of 
Ur < 24 mag arcsec"^. In their primary sample, 
IGladders et al.l detect a total of 5 arc candidates, four 
of which have length-to-width ratios r > 10. Oddly, all 
five of the arc-candidate clusters in their primary sample 
have high redshifts, z > 0.64. From the EMSS statistics, 
we would expect on average ~ 2 arc-producing clusters 
in the range 0.15 < z < 0.6 over a 90 deg^ field, so the 
lack of such clusters in the RCS sample is unusual but 
not sig nificantlv discrep ant. 

Since IGladders et al.1 impose a surface brightness cut 
on their arc candidates, to compare with their statistics 
we need to know the source galaxy numbei' density as 
a function of surface brightness. iFontana et al.l |2000) 
do not give surface brightnesses for the galaxies in their 
Hubble Deep Field catalog, but they do provide mag- 
nitudes and half-light radii for galaxies observed in the 
NTT Deep Field. From their catalog, we derive number 
densities of galaxies with fin < 24 mag arcsec"^ of 6.9, 
1.7, and 9.0 ar cmin"^ for our three source redshift bins. 
IFontana et al.l do not give estimates for their incomplete- 
ness at these surface brightnesses. From figure^! we can 
see that GIF simulation lenses at z > 0.6 primarily pro- 
duce arcs from sources at ^ 2, so for simplicity we will 
only count these sources in our statistics. 

From the GIF simulation outputs at z = 0.62 and 
z = 0.82, we would expect ~ 1 arc in the RCS field 
with length-to-width ratios exceeding 10 , and distances 
from the cluster center exceeding 10". 'Gla dders et alJ 
find four such arcs around three clusters in their primary 
sample. In addition, they also report another giant arc 
at smaller radii. Dropping the cut on > 10" increases 
the expected number of arcs with length-to-width ratio 
r > 10 to ~ 2, while painting galaxies with a ~ 350 km/s 
onto the clusters increases the expected number of arcs 
to ~ 5. 
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Fig. 11. — Number weighted cross-section at redshift z = 0.82 
for sources at Zs = 2. The strong rise out to the largest masses 
may indicate that the dominant lenses are not represented in the 
GIF volume at high redshifts. 



In addition to the arc candidates included in their pri- 
mary sample, IGladders et al.1 lj2003j) also report the dis- 
covery of three arc-producing clusters found in a follow- 
up survey of clusters with redshifts z > 0.95. Two of the 
clusters have arcs at small radii, however the third, RCS 
2319.9-^0038, has multiple arcs at > 10". 

The RCS statistics may indicate an excess of giant 
arcs observed at high redshifts, in contrast to the low- 
redshift agreement found between EMSS and the ex- 
pected GIF arc statistics. The excess appears to be 
due to the presence of extremely massive and/or concen- 
trated clusters in the RCS fields which are not present 
in the GIF volume. For example, RCS 2319.9-1-0038 has 
^ lO^^-^/i-^Mg enclosed within the central 20", while 
the most massive cluster in the z = 1 GIF volume has 
less than twice that mass inside its virial radius! The 
high fraction of multiple-arc systems in the RCS sam- 
ple further supports the notion argued bv lGladders et al.1 
that a few "super-lenses" are responsible for much of the 
optical depth, although a full treatment of multiple arc 
systems must take into consideration source clustering, 
which is cer tainly important in at least some multi-arc 
clusters fe.g lSand et"al]l2002|) . 

Unfortunately, these rare clusters are precisely the 
clusters for which the small volume (~ lO^ft.'^Mpc'^) 
we use is deficient in describing the much larger volume 
(> 10^/i~^Mpc'^ comoving) probed by the RCS sample. 
Massive objects like RCS 2319.9-H0038, which evidently 
occur only once or so per RCS volume, will almost cer- 
tainly never appear in a volume two orders of magnitude 
smaller. Indeed, examination of the giant arc optical 
depth as a function of mass at the relevant redshifts (Fig- 
ure indicates that we are missing the objects which 
dominate the optical depth. Accordingly, we cannot say 
at present whether the RCS results are consistent or not 
with the ACDM model. 

4. COMPARISON WITH PREVIOUS WORK 

We find rough agreement with the optical depths re- 
ported by B98 for giant arcs with length-to-width ratios 



r > 10. Nevertheless, we do not find that this opti- 
cal depth leads to an order of magnitude deficit in the 
number of arcs expected in EMSS clusters compared to 
the observed number. Part of the difference between our 
conclusions and those of B98 stems from the difference in 
estimates for cluster number counts. We have estimated 
this using d i rectly the EMSS sky coverage reported in 
iHenrv et all l)1992|) . while B98 assume a number density 
for EMSS clusters corresponding to about twice the sky 
density we have estimated. In addition, B98's estimate 
for the source density is lower than our estimate (based 
on counts from the Hubble Deep Fields) also by a factor 
of two. This combined with our somewhat larger opti- 
cal depths accounts for the difference in our conclusions. 
One puzzling difference between our results is that B98 
find a much steeper dependence of optical depth on arc 
length-to-width ratio. We have quoted statistics for arcs 
with r > 10. For lower thresholds, e.g. r > 7.5, our 
optical depths increase by a factor ~ 2. In contrast, B98 
find that the optical depth for r > 7.5 is nearly an order 
of magnitude larger than that for length-to-width ratio 
of 10. Theoretically, one might expect the cross section 
to scale like a oc lc5{ii), where Ic is the caustic length 
and 5 is the distance from the caustic for which sources 
are magnified by Near folds, the magnification scales 
like ^ oc while ne ar cusps the magnification scales 

like /i oc (5~^ (Schneid er et al.l ll9921. Naively, assum- 
ing r is proportional to near critical lines, one might 
expect the cross section to scale like a oc l/r^~^, consis- 
tent with the moderate rise in r with declini ng r which 
we fi nd. We al so note that o ther workers ( Qguri et alJ 
120031 (Wambsg anss et al.l2003y) find r dependence similar 
to ours in their simulations, and certainly the observed 
counts of giant arcs do not show such a steep dependence 
on length-to-width ratio. 

We do not reproduce t he enhanced optical depths of 
iWambsganss et aP 1)20031) . and in particular we do not 
find as steep a redshift dependence. For example, their 
reported optical depth at Zs = 1.5 of r = 2.3 x 10~^ is 
a factor of 3 larger than our value. It is possible that 
part of this discrepancy could be attributed to their use 
of arc magnification as a proxy for length-to- width ratio; 
as shown in figure [21 this could overestimate the cross 
section by a factor of 2-3. 

The weaker dependence on source redshift is more diffi- 
cult to explain. The main difference between our strate- 
gies was that we only included material within ~ 100 
Mpc of the cluster and did not integrate along the entire 
light cone. However, given the expected fluctuations of 
only ~ 1% of the critical surface density, it seems un- 
likely that large scale structure could affect the lensing 
cross-section significantly. It is possible to significantly 
enhance the cross-section with small fluctuations in sur- 
face density for lenses with Einstein radii much smaller 
than the scale radius, where the surface density profile 
is extremely shallow, but this is not the regime in which 
we expect to see giant cluster arcs. We have explicitly 
checked that the addition of extra surface density at the 
level of 5k ~ 3% changes the giant arc cross section by 
< 5% for the massive clusters which dominate the optical 
depth. 

The simulations of IWambsganss et aD had signifi- 
cantly higher force resolution than the GIF simulations: 



9 



iWambseanss et all used N-body simulations with ~ 3 
kpc force resolution, while the GIF simulations had force 
resolution of 20 kpc. However, 20 kpc is still well 
within the typical radii at which these giant arcs form, 
and we have argued in ij2 2. 51 that this resolution is suffi- 
cient to calculate the total optical depth to forming giant 
arcs. At present, we do not have a good understanding 
for the different behavior with redshift but our investi- 
gations of the effects of central galaxies showed that it is 
possible that the redshift behavior is sensitive to which 
arcs are considered to be giant cluster arcs. 

The enhanced cross-sections found by iQguri et al.l 
|)200.1D are also puzzling, since the main effect seems to 
come from considering steeper profiles in the inner re- 
gion (~ 10 — 20 kpc), well inside typical Einstein radii 
(~ 50 — 100 kpc). Again, this could plausibly be a sensi- 
tivity in the total cross-section to the details of the cen- 
tral region that is in reality not relevant to the problem 
o f giant cluster ar c s. 

iGladders et al.l l)2003D have suggested that "super- 
lenses", objects with low mass but large cross section, 
may be responsible for much of the optical depth. Their 
argument for low masses hinges on the assumption that 
the cross-section for giant arcs scales linearly with the 
mass, but fig.0]shows that a more reasonable estimate is 
a cx M^~'^. With this steeper scaling it is plausible that 
"superlenses" are primarily simply more massive clus- 
ters. Weak lensing mass maps at larger radii or Sunyaev- 
Zeldovich measurements can answer the question of 
whether high-redshift lenses like RCS 2319.9-^0038 are 
supermassive (e.g. M ~ 10^^h~^MQ) or not. If these 
objects have lower masses (e.g. fewx 1 0^^/i~^Mfn), then 
their la rge Einstein radi i are p uzzling. iMeneghetti et all 
l)2003c|) and lTorri et al.l l|20fl3f) have suggested that clus- 
ter mergers may account for a significant fraction of 
the giant arc optical depth. However, note that the 
critical lines tr aced out by the multiple-arc systems of 
IGladders et alJ appear relatively round, in contrast to the 
highly elongated critical lines expected of low mass clus- 
ters whose giant arc cross sections have been enhanced 
by on going mergers ( Mcneghctti ct al. 2003c; Torri ct all 
l2003j) . Another possibility is that "super-lenses" may be 
ordinary mass, but super-concentrated. An example of 
this may be CL0024-I-1654, for wh ich a recent combi ned 
strong- and weak-lensing analysis l)Kneib et alJl2003|) in- 
dicates a concentration C200 ~ 22. Such high concen- 
trations are extrem ely anomalous for s imulated CDM 
galaxy clusters (e.s iBullock et al.ll200lh . If the RCS 
super-lenses turn out to have similarly deviant concentra- 
tions, and such objects are commonplace, their existence 
could pose a serious challenge to our understanding of 

http:/ /www.cfht.hawaii.odu/Science/CFHLS 
^ http://www.sdss.org 



structure formation in the ACDM model. 

5. DISCUSSION AND CONCLUSIONS 

At present, there appears to be no discrepancy between 
observed giant arc statistics and predictions from the 
ACDM mod el at low redshif ts z < .6. Imaging of EMSS 
clusters by iLe Fevre et al.l l|1994D and Luppino et alJ 
il999*) indicates that ~ 900 giant arcs are expected over 
the full sky from X-ray bright clusters at these redshifts, 
while ray-tracing through clusters taken from the GIF 
simulations leads us to expect ^ 1000 giant arcs. At 
higher redshifts z > 0.6, the observed abundance of close- 
in arcs with d < 10" observed by the LCDCS and RCS 
appears consistent with expectations, when the effects of 
the central galaxies are included in the gravitational po- 
tential. However, there may be an excess of wide separa- 
tion {d > 10") giant arcs at high redshift. This putative 
excess appears to be due to extremely massive or unusu- 
ally concentrated clusters not present in the (admittedly 
limited) simulation volumes we have employed. 

Looking for ward, upcomin g X-ray cluster surveys like 
MACS (Ebeli ng et all 120011 and wide -area surveys like 
the CFHT Legacy Survey'' and the Sloan Digital Sky Sur- 
vey^ can be expected to improve the statistics of giant 
arcs on the sky. For example, the RCS-2 survey covers an 
area of 830 deg^ and is expected to produce ~ 50 — 100 
new arcs (M. Gladders, private communication). Deep 
surveys can also begin to test CDM predictions for giant 
arc properties. As we have discussed earlier, CDM clus- 
ters tend to produce arcs along (but aligned orthogonal 
to) their major axes. This is modified somewhat by the 
presence of central cD galaxies in the clusters, however 
for reasonable velocity dispersions we expect the arc an- 
gular distribution to be anisotropic. Multi-arc systems 
may provide a probe of halo triaxiality at small radii, 
complementary to we ak lensing studies at large radii (e.g. 
iHoekstra et~al]l2003|) . 
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APPENDIX 
ON LIPS AND NAKED CUSPS 

In section H2 2.41 the addition of a centrally concentrated mass was shown to have a significant effect on where giant 
arcs can form. We can understand this behavior as follows. On the critical curves, the magnification matrix becomes 
singular, causing sources to become highly magnified along one direction. The highly elongated direction, however, is 
not in general parallel to the critical line itself. In some instances, the magnified direction can be perpendicular to 
the critical line, i.e. it can point along V/i. Such cases clearly are not conducive to the formation of giant arcs, as 
the magnification diminishes so quickly along their lengths that these arcs are unable to acquire 'giant' dimensions. 
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Instead, giant arcs form preferentially at locations where the magnified direction is aligned with the tangent to the 
critical curve. 

The direction along which arcs are elongated is determined by the orientation of the local shear. For our purposes, 
the shear has two principal sources: the tidal gravity associated with the radial density gradient, and the triaxiality 
and substructure of the dark matter distribution. The former produces shear which is azimuthally oriented with 
respect to the central mass concentration, while the latter lead to shear aligned with the principal axes of the density 
distribution. The competition between these two effects determines whether arcs are primarily tangential or 'Cartesian', 
which determines whether the caustic structure is that of lips, or radial and tangential caustics. 

N-body simulations in the cold dark matter model tend to produce clusters which have shallow radial profiles and 
strong triaxiality. Average radial profiles exhibit density cusps fe.g lNavarro et al]ll997j) which, in projection, give 
surface densities rising only logarithmic ally at small r adii. Typ ical clusters are also quite triaxial, with axis ratios on 
average near 0.5 with large scatter fe.g lJing fc Sutoi r2002: Th omas et al.lll998j) . Because of this, a large fraction of 
GIF clusters have lips caustics, although this depends in part on projection. Clusters viewed along their major axes 
have larger central surface densities and less ellipticity than other projections, meaning that clusters can have radial 
and tange ntial c austics viewed along certain projections, and lips when viewed in other orientations. 

In figure fST2l we plot the caustics and critical lines for an elliptical NFW profile, and lensed images of representative 
sources. The top panel illustrates typical caustic patterns for the pure DM simulations. There are two critical lines, 
pointing along and orthogonal to the cluster major axis. However, as the lensed arcs indicate, the magnified direction 
along the outer critical line tends to point up and down. Because of this, the corresponding caustics in the source 
plane take the form of lips catastrophes, with two cusps each, corresponding to the two locations on each critical line 
where the highly magnified direction becomes tangent to the curve ( S chneider, et al. 1992). This is why long, thin arcs 
preferentially form at the ends of the critical lines in this case. 

The caustic structure changes as the central mass concentration increases, shown in the bottom two panels of the 
figure. As the central mass is increased, the tangential shear grows commensurate to the ellipticity-induced shear, 
and the two critical lines approach each other. When they meet, their corresponding caustics meet at a hyperbolic 
umbilic catastrophe (Schneider et al. 1992), shown in the middle panel of figure IXT2I Near hyperbolic umbilics, the 
source becomes highly magnified in both directions, which again is not favorable for producing long, thin arcs. For 
higher central concentrations, the critical curves become tangential and radial; that is, their magnified directions tend 
to point in the azimuthal and radial directions, for the outer and inner curves, respectively. The near alignment of the 
magnified direction with the curve tangent, for t he ta ngential critical line, allows long thin arcs to form all along its 
length, as depicted in the bottom panel of figure IXT2I 
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Fig. A12. — Critical lines and caustic structure as a function of central mass concentration. The top panel depicts lensing by an elliptical 
NFW profile, while the bottom two panels show the effect of adding increas ingly massive c entral galaxies, modeled as Jaffe profiles. Critical 
lines, caustics, and arcs were computed in this figure using the software of lKeetonI 120011) . 
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